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Abstract 

Herein we dennonstrate how nanojunctions between lysosonnes and 
sarcoplasnnic reticulunn (L-SR junctions) serve to couple lysosonnal activation to 
regenerative, ryanodine receptor-nnediated cellular Ca^+ waves. In pulnnonary 
artery snnooth nnuscle cells (PASMCs) it has been proposed that nicotinic acid 
adenine dinucleotide phosphate (NAADP) triggers increases in cytoplasnnic Ca 
^+ via L-SR junctions, in a nnanner that requires initial Ca^+ release fronn 
lysosonnes and subsequent Ca^+-induced Ca^+ release (CICR) via ryanodine 
receptor (RyR) subtype 3 on the SR nnennbrane proxinnal to lysosonnes. L-SR 
junction nnennbrane separation has been estinnated to be < 400 nnn and thus 
beyond the resolution of light nnicroscopy, which has restricted detailed 
investigations of the junctional coupling process. The present study utilizes 
standard and tonnographic transnnission electron nnicroscopy to provide a 
thorough ultrastructural characterization of the L-SR junctions in PASMCs. We 
show that L-SR nanojunctions are pronninent features within these cells and 
estinnate that the junctional nnennbrane separation and extension are about 15 
nnn and 300 nnn, respectively. Furthernnore, we develop a quantitative nnodel of 
the L-SR junction using these nneasurennents, prior kinetic and specific Ca^+ 
signal infornnation as input data. Sinnulations of NAADP-dependent junctional 
Ca^+ transients dennonstrate that the nnagnitude of these signals can breach 
the threshold for CICR via RyR3. By correlation analysis of live cell Ca^+ 
signals and sinnulated Ca^+ transients within L-SR junctions, we estinnate that 
"trigger zones" connprising 60-1 00 junctions are required to confer a signal of 
sinnilar nnagnitude. This is connpatible with the 1 10 lysosonnes/cell estinnated 
fronn our ultrastructural observations. Most innportantly, our nnodel shows that 
increasing the L-SR junctional width above 50 nnn lowers the nnagnitude of 
junctional [Ca^+] such that there is a failure to breach the threshold for CICR via 
RyR3. L-SR junctions are therefore a pre-requisite for efficient Ca^+signal 
coupling and nnay contribute to cellular function in health and disease. 
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Introduction 

The importance of cytoplasmic nanoj unctions to cellular signaling 
and thus to the modulation of cell function was recognised several 
decades ago^ -, henceforth the extent to which cellular nanospaces 
may contribute to the regulation of cell function received little 
attention. Nevertheless there is now a growing recognition of the 
widespread occurrence and functional significance of cytoplasmic 
nanospaces in cells from species across several kingdoms^"'^. 

In this respect, membrane-membrane junctions between lysosomes 
and the sarcoplasmic reticulum (L-SR junctions) are of particular 
interest; not least because they have been hypothesized to couple 
calcium signaling between these two organelles^^'^^ 

That L-SR junctions may play an important role in cellular Ca^^ 
signaling was uncovered through early studies on the Ca^+ mobi- 
lizing messenger nicotinic acid adenine dinucleotide phosphate 
(NAADP)^% which demonstrated that NAADP released Ca^+ from 
a store other than the sarco/endoplasmic reticulum (S/ER)^\ that 
could then trigger further Ca^^ release from the S/ER by Ca^^- 
induced Ca^^ release (CICR)^^"'^. A major advance in our under- 
standing was then provided by the demonstration that the NAADP- 
released Ca^+ was from an acidic lysosome-related store"^'^^'^'^ in a 
manner that requires two pore segment channel subtype 2 (TPC2)-^\ 
However, studies on pulmonary arterial smooth muscle cells 
(PASMCs) had also identified a significant specialization, namely 
L-SR nanoj unctions. It was hypothesized not only that these 
nanojunctions were necessary for coupling between lysosomes 
and the SR, but that they could both coordinate and restrict their 
relationship to the SR by preferentially targeting ryanodine recep- 
tors while excluding inositol 1,4,5-trisphosphate (IP3) receptor s^^^'^^ 
Importantly, NAADP-dependent Ca^^ bursts primarily arise from 
lysosomes in the perinuclear region of PASMCs and appear to pro- 
mote rapid, local Ca^+ transients that are of sufficient size to activate 
clusters of SR resident ryanodine receptor subtype 3 (RyR3) and 
thus initiate, in an all-or-none manner, a propagating global Ca^+ 



The specialization of the proposed L-SR junction is likely impor- 
tant in smooth muscle cell physiology, e.g., in vasomotion, and 
in the recycling of organelles and programmed cell death by 
autophagy-^ However, L-SR junctions may also make as yet unfore- 
seen contributions to vascular pathologies as highlighted by the 
fact that Niemann-Pick disease type CI results, in part, from dys- 
regulation of lysosomal Ca^^ handling-- and is known to precipitate 
portal hypertension-^ while other lysosomal storage diseases (e.g., 
Pompe and Gaucher disease) accelerate pulmonary arterial hyper- 
tension-'^'-^ Moreover, observed hypertension is often associated 
with dysfunction of cholesterol trafficking-^', increased plasma cho- 
lesterol levels, vascular lesion formation, atherosclerosis/thrombo- 
sis and medial degradation-^ -^. Therefore lysosomal Ca^^ signaling 
is of considerable clinical interest. That L-SR junctions may be of 
further significance to pathology is also evident, for example, from 
the fact that in the pulmonary artery smooth muscle L-SR junctions 
underpin Ca^^ waves initiated by endothelin 1, the levels of which 
are elevated in pulmonary hypertension, systemic hypertension and 
atherosclerosis-^'-'^. An understanding of how specific Ca^^ signals 



are functionally initiated therefore has important translational 
implications. 

Lysosomal Ca^+ regulation has been of current interest in several 
recent research and review articles (e.g.,^""^^), yet a quantitative elu- 
cidation of the mechanism by which Ca^^ signals are generated by 
the endolysosomal system is still lacking. Given the likely impor- 
tance of L-SR junctions to Ca^+ signaling in health and disease, we 
sought to determine whether this nano-environment would indeed 
be able to effectively generate the previously observed NAADP- 
induced Ca^^ signals. In particular, we set out to verify the following 
hypotheses in PASMCs: (1) L-SR nanojunctions should be observ- 
able in the ultrastructure of these cells, (2) NAADP induces suf- 
ficient Ca^+ release from the lysosome to initiate activation of RyR3 
embedded in the junctional SR, and (3) the combined effect of acti- 
vation of L-SR nanojunctions in a cytoplasmic "trigger zone" deter- 
mines the threshold of global [Ca^^]. for the biphasic release process. 

Due to the minute spatial scale of the nanojunctions generating the 
primary Ca^^ signals, accurate investigation of dynamic signaling 
within these spaces cannot be achieved with currently available 
instrumentation. To overcome this issue, we took an integrative 
approach by combining our own electron microscopy of L-SR 
nanojunctions, existing kinetic data on the Ca^^ transporters and 
buffers, and prior knowledge of the NAADP-induced Ca^^ signal 
features into a quantitative model of a typical L-SR nanojunction. 
This model is based on stochastic simulations of intracellular Ca^^ 
diffusion by Brownian motion implemented using the particle sim- 
ulator MCelp5-^l 

In the present manuscript we have verified the existence of L-SR 
nanojunctions within the ultrastructure of PASMCs, and shown that 
lysosomes can release sufficient Ca^^ to activate CICR via RyR3 
clusters embedded in the junctional SR. Perhaps most importantly, 
we show that L-SR coupling is determined both by the integrity of 
L-SR junction on the nanoscale and the quantal summation of Ca^+ 
release from multiple, activated junctional complexes. 

Materials and methods 

All the experiments and procedures were carried out in accordance 
with the guidelines of the University of British Columbia Animal 
Care Committee, and in accordance with the United Kingdom 
Animals (Scientific Procedures) Act 1986. The Research Ethics 
Committees from both these institutions specifically and indepen- 
dently approved this study. 

Electron microscopy 

An adult male Wistar rat (about 12 months old, weighing 220-250 g, 
obtained from Charles River and housed in the institutional ani- 
mal facility at the University of British Columbia Child and Family 
Research Institute under standard animal room conditions of 12 h 
light- 12 h dark cycles, at 25 °C, with two animals in a cage) was 
anesthetized with 3 ml of Thiotal (Vetoquinol, Lavaltrie, QC, Can- 
ada). For the purpose of this phase of the study, one representative 
healthy animal usually suffices. This choice also minimizes unnec- 
essary animal sacrifice. The thoracic cavity was opened and flooded 
with warm physiological saline solution (PSS) containing 145 mM 
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NaCl, 4 mM KCl, 1 mM MgCl^, 10 mM HEPES, 0.05 mM CaCl^, 
and 10 mM glucose (pH 7.4; all chemicals were purchased from 
Sigma- Aldrich Canada unless otherwise stated). After extraction of 
the heart and lungs and their placement in HEPES buffer, several 
rings from the primary and secondary branches of the pulmonary 
artery were dissected and immediately immersed in fixative solu- 
tion. The primary fixative solution contained 2.5% glutaraldehyde 
(Ted Pella, Redding, CA, USA) in 0.1 M sodium cacodylate buffer 
(Canemco & Marivac, Gore, QC, Canada) at room temperature. 
The artery rings were then washed three times in 0.1 M sodium 
cacodylate (30 min in total). In the process of secondary fixation, 
the tissue rings were fixed with 1% OsO^ (Ted Pella, Redding, CA, 
USA) in 0.1 M sodium cacodylate buffer for 1 h followed by three 
10-minute washes with distilled water and en bloc staining with 2% 
uranyl acetate (Ted Pella, Redding, CA, USA. The samples were 
then dehydrated in increasing concentrations of ethanol (25, 50, 75, 
80, 90, and 95%). In the final process of dehydration, the samples 
underwent 3 washes in 100% ethanol. The artery rings were then 
resin-infiltrated in increasing concentrations (30, 50, and 75% in 
ethanol) of a 1:1 mix of Epon and Spurr's resins (Canemco & Marivac, 
Gore, QC, Canada). The infiltration process was completed by three 
passages in 100% resin. All of the ethanol dehydration and resin 
infiltration steps were carried out by using a laboratory microwave 
oven (Pelco 3441, Ted Pella, Redding, CA, USA). The blocks were 
finally resin-embedded in molds and polymerized overnight in an 
oven at 60°C. 

2D imaging 

For standard (2D) electron microscopy imaging, 80-nm sections 
were cut from the embedded sample blocks on a Reichert Ultracut-E 
microtome (Leica Microsystems, Vienna) using a diamond knife 
(Diatome, Biel, Switzerland) and were collected on uncoated 100- 
and 200- mesh copper grids (hexagonal or square meshes; Ted 
Pella, Redding, CA, USA). The sections were post-stained with 1% 
uranyl acetate (Canemco & Marivac, Gore, QC, Canada; pH not 
recorded) and Reynolds lead citrate (Fisher Scientific Company, 
USA; pH not recorded) for 12 and 6 minutes, respectively. Electron 
micrographs at various magnifications were obtained with a Hitachi 
7600 transmission electron microscope (Hitachi High Technolo- 
gies, Japan) at 80 kV. 

Lysosomes in these images were identified according to their well 
established appearance and features: they are bound by a single 
lipid-bilayer membrane, with a granular, more or less uniform lumi- 
nal matrix that is more electron dense than the surrounding cyto- 
sol. Secondary lysosomes may also contain less granular structures 
within the finer matrix. Moreover, lysosomes are normally dis- 
tinguishable from endosomes by their larger size, hence we set a 
threshold "diameter" of > 200 nm for acceptance of a lysosome, 
below which all vesicles were excluded. 

Tomography (3D imaging) 

To obtain electron microscopic tomograms, we cut 200-nm-thick 
sections from the same sample blocks used for standard imaging. 
The sections were then collected on Formvar (Canemco & Marivac, 
Gore, QC, Canada) coated slot copper grids (Ted Pella, Redding, 
CA, USA) and post-stained with 1% uranyl acetate and Reynolds 
lead citrate for 20 and 10 minutes, respectively. We surveyed the 



sample sections using a FEI Tecnai G2 200 kV transmission elec- 
tron microscope (FEI Technologies, Eindhoven, The Netherlands) 
and identified regions of interest containing L-SR junctions. We 
then acquired tomograms of several of those regions by taking 2D 
scans through the sample every 5° of inclination as it was tilted 
between -60° and +60° with respect to horizontal. The scans were 
reconstructed with the Inspect3D (FEI Technologies, Eindhoven, 
The Netherlands) software tools and structures of interest, primar- 
ily lysosomal and SR membranes in the same cellular neighbour- 
hood, were traced out using Amira (FEI Technologies, Eindhoven, 
The Netherlands) software. 

Image and data analysis 

The images of the samples were further processed using GIMP 
(GNU Imaging Manipulation Program, open source, available at 
gimp.org) to enhance membrane contrast in the interest of improv- 
ing our characterization of the L-SR junctions. 

The SR and lysosomal membranes were outHned, highlighted and 
measured in pixels using the Inkscape software program (open 
source, available at http://www.inkscape.org), converting the pixel 
measurements to nm using the scale bar in the recorded micro- 
graphs. By modifying the Inkscape script for measuring lengths, 
we were able to output the measurements directly into a text file, 
which we used to produce the histograms in Figure 3B-D. We used 
the Gnuplot program (open source, available at http://www.gnuplot. 
info) to produce the histograms and plots presented herein. 

Ca^+ imaging 

To obtain sufficient cell isolation to appropriately test significance, 
1-3 adult male Wistar rats (weighing 150-300 g, obtained from 
Charles River and housed in the institutional animal facility at the 
Centre for Integrative Physiology of the University of Edinburgh 
under standard animal room conditions of 12 h light- 12 h dark 
cycles, at 25 °C, with 2-4 animals in a cage) were sacrificed by cer- 
vical dislocation. The heart and lungs were removed en bloc and 
placed in physiological salt solution (PSS) of the following compo- 
sition (mM): 130 NaCl, 5.2 KCl, 1 MgCl^, 1.7 CaCl^, 10 glucose, 
and 10 HEPES, pH 7.45 (all chemicals were purchased from Sigma- 
Aldrich). Single arterial smooth muscle cells were isolated from 
second-order branches of the pulmonary artery. Briefly, arteries 
were dissected out and placed in low Ca^+ solution of the following 
composition (mM): 124 NaCl, 5 KCl, 1 MgCl^, 0.5 NaH^PO^, 0.5 
KH^PO, 15 NaHC03, 0.16 CaCl^, 0.5 EDTA, 10 glucose, 10 taurine 
and 10 HEPES, pH 7.4. After 10 min the arteries were placed in 
the same solution containing 0.5 mg/ml papain (Ruka) and 1 mg/ml 
bovine serum albumin (Sigma) and kept at 4°C overnight. The fol- 
lowing day 0.2 mM 1,4-dithio-DL-threitol (Sigma) was added to 
the solution, to activate the protease, and the preparation was incu- 
bated for 1 h at room temperature (22°C). The tissue was then 
washed at 3°C in fresh low Ca^^ solution without enzymes, and 
single smooth muscle cells were isolated by gentle trituration with 
a fire-polished Pasteur pipette. Cells were stored in suspension at 
4°C until required. 

PASMCs were incubated for 30 min with 5 ^M Fura-2-AM in Ca^^- 
free PSS in an experimental chamber on a Leica DMIRBE inverted 
microscope and then superfused with Fura-2 free PSS for at least 
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30 min prior to experimentation. Intracellular Ca^^ concentration 
was reported by Fura-2 fluorescence ratio (F340/F380 excitation; 
emission 510 nm). Emitted fluorescence was recorded at 22°C with 
a sampling frequency of 0.5 Hz, using a Hamamatsu 4880 CCD 
camera via a Zeiss Fluar 40x, 1.3 n.a. oil immersion lens and Leica 
DMIRBE microscope. Background subtraction was performed 
on-line. Analysis was done via Openlab imaging software (Impro- 
vision, UK). 

NAADP was applied intracellularly in the whole-cell configuration 
of the patch-clamp technique, and in current clamp mode (/ = 0) 
as described previously The pipette solution contained (in mM): 
140 KCl, 10 HEPES, 1 MgCl^ and 5 /jM Fura-2, pH 7.4. The seal 
resistance, as measured using an Axopatch 200B amplifier (Axon 
Instruments, Foster City, CA), was > 3 GQ throughout each experi- 
ment. Series resistance and pipette resistance were < 10 MQ and 
< 3 MQ, respectively. All experiments were carried out at room 
temperature (» 22°C). 

Quantitative modeling 

The main stages of the quantitative modeling approach are: 

1. The design of 3D software mesh objects (nets of intercon- 
nected triangles by which surfaces can be represented in com- 
puter graphics) representing a typical L-SR region, including 
a whole lysosome and a portion of neighbouring SR network. 
These objects are built to-scale following the ultrastructural 
characterization of the L-SR junctional regions as it results 
from our electron microscopy image analysis; this phase was 
carried out using the "3D content creation suite" Blender (open 
source, available at blender.org); 

2. The positioning of the relevant transporters on the reconstruct- 
ed membranes (TPC2 complexes on the lysosome, SERCA2a 
and RyR3 on the SR) according to information gathered from 
the literature on their typical membrane densities and the im- 
plementation of the transporter's known kinetics and multi- 
state models and of the ion diffusivities (Ca^^ and mobile Ca^^ 
buffers); 

3. The simulation of molecular Brownian motion in the cytosol 
by random walk algorithms; this phase was performed by writ- 
ing appropriate code for the stochastic particle simulator MCell 
(freely available at http://www.mcell.org)^^"^^. In a nutshell, 
MCell reproduces the randomness of the molecular trajecto- 
ries, of the ion transporter flickering and of the relevant chemical 
reactions by probabilistic algorithms governed by pseudo- 
random number generators (iterative mathematical algorithms, 
which produce a pseudo-random sequence of numbers once 
initiated by a given number called seed). This enables the sim- 
ulation of a number of microphysiological processes, all sto- 
chastically different from one another. The average outcome 
of the processes, to mimic the instrumental output during 
experimental measurements, is obtained by taking the average 
of a desired quantity, e.g., [Ca^+], over a large number of simu- 
lations all initiated with a different seed; 



4. The measurement of simulated [Ca^^] in the L-SR junctions 
from the process of Ca^+ release via TPC2-related signaling 
complex on the lysosome and SR Ca^+ uptake by the SERCA2a 
pumps, and static as well as dynamic visualization of the simu- 
lations; this stage is part of MCell's output. 

Results 

The work presented in this article is based on earlier Ca^^ fluores- 
cence imaging observations by optical microscopy, of which the 
data in Figure 1 is a representative example, and immunofluores- 
cence results. These studies led to the proposal that, for the lysoso- 
mal Ca^^ release to trigger CICR, L-SR nanojunctions are required 
and that they consist of apposing patches of lysosomal and SR 
membranes separated by a narrow space of nano- scale dimen- 
§iQj^io,i7 jj^g relevant background findings that stimulated the devel- 
opment of the work exposed here were first reported by Evans' 
group"^'^^ and are summarized in Figure 1. 

The example record in Figure 1 highlights the fact that NAADP 
appears to activate a two-phase Ca^^ signal, which can exhibit an 
identifiable "shoulder" during the initial rising phase of the signal 
(Figure 1, time point 1), followed by a second faster phase of sig- 
nal amplification (Figure 1, time point 2). Previous studies have 
demonstrated that this two-phase response results from initiation by 
NAADP of Ca^^ bursts from lysosome-related stores"^ in a manner 
that requires TPC2^^ and that Ca^^ bursts are subsequently amplified 
by CICR from the SR via clusters of RyR3^^'^^. The studies cited 
above provided calibrated estimates of the changes in intracellular 
[Ca^^] input into the model below. We propose herein that L-SR 
junctions do indeed represent cellular nanojunctions and that they 



/yy y y y| 
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WC 10nM NAADP 



Figure 1. Upper panel shows a series of pseudocolour images of 
tlie Fura-2 fluorescence ratio (F340/F380) recorded in an isolated 
pulmonary artery smooth muscle cell during intracellular dialysis of 
10 nM NAADP. Note the spatially localized 'Ca^^ burst' (time point 1). 
Lower panel shows the record of the Fura-2 fluorescence ratio 
against time corresponding to the upper panel of pseudocolours 
images; note the discrete shoulder in the rising phase of the F340/ 
F380 ratio that corresponds to the initial 'Ca^^ burst'. Scale bar: 1 0 /^m. 
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might play a role of accentuating Ca^+ gradients, akin to that of 
plasma membrane (PM)-SR junctions that are pivotal in the pro- 
cess of SR Ca^+ refilling during asynchronous [Ca^+] waves^'^. We 
furthermore hypothesize that in order for these nanojunctions to 
appropriately regulate Ca^^ signaling, they must be separated by a 
distance of approximately 20 nm and be of the order of a few hun- 
dred nm in lateral dimensions, as inferred from previous studies on 
PM-SRjunctions^O'i^ 

Ultrastructural characterization of L-SR nanojunctions 

To identify lysosomes, SR regions and L-SR nanojunctions, we 
recorded and surveyed 74 electron micrographs of rat pulmonary 
arterial smooth muscle taken from samples prepared as described 
in the Materials and methods section. The images in Figure 2 



provide a set of examples. Since we were aiming to detect L-SR 
junctions, namely close appositions of the lysosomal and SR 
membranes, immuno-gold labeling of lysosomes was prohibited, 
since this technique compromises membranes definition by electron 
microscopy to the extent that we would be unable to assess junc- 
tional architecture. Instead, in images like those in Figure 2, lyso- 
some identification was accomplished by relying on the knowl- 
edge of lysosomal ultrastructural features, which has accumulated 
over the past 50 years since the discovery of the lysosomes^^\ In 
standard (2D) transmission electron microscopy (TEM) images, lys- 
osomes typically appear as elliptical structures bound by a single 
lipid bilayer, a feature that distinguishes them from mitochondria. 
Depending on the lysosomal system stage, they also tend to have a 
more or less uniformly electron-dense interior as compared to the 




Figure 2. Representative electron micrographs of rat pulmonary artery SMC regions containing lysosomes (L), several SR cisterns, and 
including several examples of L-SR junctions (arrows). Also indicated are nuclei (N), Golgi apparatus (G), mitochondria (M), a multivesicular 
body (MVB) and extra-cellular space (ECS). Scale bars: 500 nm. Magnifications: A,C 80,000x, B, 60,000x, D, 70,000x. 
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surrounding cytosol^^ They can be distinguished from endosomes 
by their larger size and darker lumen and they differ from peroxi- 
somes, since the latter usually display a geometrically distinct and 
markedly darker structure called "crystalloid" in their luminal area. 
Moreover, it would appear that peroxisomes are found far more fre- 
quently in liver, kidney, bronchioles and odontoblasts than in other 
cell types^^^'^- (see, for example, http://www.uni-mainz.de/FB/Medi- 
zin/Anatomie/workshop/EM/EMPeroxisomE.html). 

Occasionally, organellar remnants are still visible inside these 
ovals, a characteristic that identifies them as multi-vesicular bod- 
ies ("MVB" in Figure 2A). As it is at times questionable whether 
MVB's are late endosomes or endosome-lysosome hybrids^"'^\ we 
have excluded organelles (3 in total) with such characteristics from 
our statistical count. 

From each of the relevant smooth muscle regions surveyed, we 
obtained high-resolution images of areas containing lysosomes and 
L-SR junctions in order to quantitatively characterize them (arrows 
in Figure 2 and Figure 3 A). Using a software graphics editor 



(inkscape.org) and the image scale bar as a calibration gauge, we 
measured the lysosome size, as the length of the major and minor 
axes of their elliptical 2D projections (in orange and grey, respec- 
tively, in Figure 3A), the L-SR widths, that is the distance between 
lysosomal and SR membranes at places where the two were about 
30 nm or closer to each other (in purple in Figure 3 A), and the 
L-SR extensions as a percentage ratio between the junctional SR 
and the lysosomal membranes (in turquoise in Figure 3 A). From 
these measurements, we extrapolated the 3D junctional SR exten- 
sion, both as a percentage of the lysosomal surface and as a length 
in nm. The histograms displayed in Figure 3B-D visually sum- 
marize the data collected from the image analysis. The mean and 
standard deviation values of the measured parameters are reported 
in Table 1 . 

Estimates of the various parameters gathered in this phase of the 
study were used as a basis to build a 3D software object to represent a 
typical L-SR junction. This reconstruction was then used to design 
the simulations mimicking Ca^^ diffusion in the L-SR nanojunc- 
tions, as is described below. 
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Figure 3. A, High magnification (150,000x) electron micrograph of a region of Figure 2B containing 3 L-SR junctions (arrows); coloured 
tracings as shown were used to measure lysosome dimensions, L-SR widths and extensions. Scale bar: 100 nm. B-D, Histograms showing 
distribution of several relevant lysosomal and L-SR junctional parameters, used to characterize the junctions and inform the quantitative 
model. B, lysosomal dimensions as major and minor axes of oval shapes in micrographs; C, L-SR junctional width; D, percentage apposition 
between junctional SR and lysosome perimeter as projected in 2D micrographs. 



Page 7 of 24 



FIOOOResearch 2014, 3:93 Last updated: 08 AUG 2014 



Table 1. L-SR junction characterization parameters. 



Parameter 


l\/lean ± SD 


Notes ^^^^^^^ 


lysosome minor axis 


(325 ± 65) nm 


n = 29 "^^^^ 


lysosome major axis 


(398 ± 91) nm 


n = 29 


L-SR width 


(16 ± 7) nm 


/I = 41 


L-SR overlap 


(23 ± 13)% 


n = 29 


L-SR lateral 
extension 


(262 ± 209) nm 


calculated from data 
above 


3D L-SR overlap 


(13 ± 15)% 


extrapolated from 2D 
measurements 



Mean and standard deviation values of L-SR junction parameters, calculated 
from data as in Figure 3B-D. 



Tomography 

To gather more direct information on the 3D morphology of L-SR 
junctions, we acquired a set of tomograms of those regions from 
the same sample blocks used to obtain the images in Figure 2. In 
Figure 4, we report snapshots from one of the tomograms; in these 
stills, we have also traced out parts of one lysosome and the closely 
apposed SR region that together form a L-SR junction. These 
tomograms are very helpful in clarifying the detailed morphology 
of L-SR junctions and informing on the spatial variability of the 
SR network. For example, while the SR segment shown in a single 
2D tomographic scan (Figure 4A) would appear to be continuous 
and part of a large SR compartment, it actually branches out into 
narrower cisterns as revealed by 3D tomographic reconstruction 
(Figure 4B). Furthermore, it is interesting to observe the fact that 
one extension of the SR appears to couple with multiple organelles. 
Thus, the 3D views generated by tomography are paramount for 
demonstrating the presence of a true junctional complex and for 
the design of a prototypical L-SR environment as a software mesh 
object, on which we may simulate the NAADP-mediated localized 
Ca^^ release. 



Quantitative model 

The model aims to verify whether NAADP-induced Ca^^ release 
from the lysosomal system could be responsible for the localized 
Ca^+ signal preceding the global Ca^^ wave (see Figure 1), which 
triggers a propagating wave by CICR via RyR3s localized at L-SR 
junctions, as predicted by earlier observations "\ 

Generation of the signal's "shoulder" by lysosomal Ca^+ 
bursts 

To understand the generation of the signal "shoulder", such as that 
observed at time point 1 in the lower panel of Figure 1, by Ca^+ 
bursts within L-SR nanojunctions, let us note that its magnitude 
corresponds to the difference in [Ca^^] between the resting level of 
» 100 nM^^ prior to NAADP stimulation (up to the point, at which 
NAADP enters the cytoplasm under the whole-cell configuration 
(WC) in Figure 1) and the value of » 400 nM^^, corresponding to 
[Ca^^]. at time point 1 in the example record shown in Figure 1, 
lower panel (these concentration values were obtained via a stand- 
ard calibration procedure^^). This leads to a A[Csl^^]^^^^^^^^ of approx- 
imately 300 nM^^. Let us now estimate the number of lysosomes 
required to generate a A[C3^^]^^^^^^^^ of such magnitude in a PASMC, 
given the following assumptions: 

1. That the luminal [Ca^+] of a lysosome, [Ca^+Jj^^, is in the range 
of 400-600 jjM, as determined in mouse macrophages^^ and 
that it is homogeneous across the lysosomal population; 

2. That the Ca^^ release rate during bursts is » 10^ ions/s (based 
on other experimenters' findings^\ but see next section for a 
detailed analysis on the rate time- variation) and that it becomes 
negligible for values of [Ca^+Jj^^ below » 80 juM, as suggested 
by the single channel kinetics of TPC2 signaling complexes in 
lipid bilayer studies^^; 

3. That lysosomes are spheroids with dimensions gathered 
from our EM characterization (see Figure 2, Figure 3, Figure 4, 




Figure 4. A, Snapshot from a TEM tomogram of a L-SR region of rat pulmonary artery smootli muscle, illustrating, among other things, a 
single SR extension apparently forming junctions with several lysosomes. Magnification: 62,000x. B, Same snapshot shown in A, but with a 
lysosome (orange) and a portion of SR (turquoise) partially traced out in 3D. These pseudocolour tracings underscore how the SR can appear 
as a large cistern in a given plane, but can actually branch out in different directions when viewed in 3D. Scale bars « 100 nm. 
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and Table 1), and hence that their volume is on average 
y=2.2x 10-^^ L; 

4. That a smooth muscle cell cytosolic volume can be calculated 
as y = 2.4 X 10"^^ L by modeling a cell as a 130-^m-long 
cylinder, 6 jjm in diameter^^\ and accounting for nuclear, SR, 
mitochondrial and lysosomal volumes as described in the sup- 
plementary material. 

With these assumptions accepted, from points 1. and 2. above we 
gather that a lysosome can release a potential A[Ca^+]jy^ = (400 to 
600) juM - 80 = 3.2 x 10 to 5.2 x 10 ^ M into'the cytosol 
(as we elaborate in the next section, this should take » 0.03 s, a time 
frame that ensures that released Ca^^ can be considered unbuffered). 
Each lysosome contribution to [Ca^^]., A[Ca^^].jy^, can then simply be 
calculated by taking the lysosome-to-cell volume ratio into account: 



A[Ca^^]„, =A[Ca^n,ys77^-3 to 5 nM (1) 

Therefore, we can calculate the number of lysosomes that may con- 
tribute to the magnitude of A[C3^^]^^^^^^^^ as 



— - 60 to 1 00 lysosomes (2) 

In summary, between 60 and 100 lysosomes would be necessary 
(and possibly sufficient) to provide a A[Csl^^]^^^^^^^^ of 300 nM, 
which is typically observed during the localized Ca^^ release phase 
of the NAADP-induced Ca^+ signal. 

How many lysosomes do we actually expect to be in a PASM cell of 
our sample tissue? We can obtain a rough estimate of this number 
from the TEM imaging we performed for this study. In each 80-nm- 
thick TEM sample section, we see 5-10 lysosomes/cell. Lysosomes 
are predominantly localized to the perinuclear region of the cyto- 
plasm, but seldom in the subplasmalemmal area, consistent with 
previous observations by optical microscopy ^^^'^^ If we simplify the 
geometry of a typical smooth muscle cell to a 130-/im-long cylinder 
with radius of 6 jjm, and considering that lysosome radii are around 
180 nm, it is reasonable to assume that separate sets of lysosomes 
would be observable in TEM images taken at distances into the sam- 
ple of about 360 nm from each other. Neglecting the slices within 
one lysosomal diameter of the cylinder surface — given the near 
total lack of observed lysosomes in those subvolumes — then our 
images suggest that we can expect a total of about 110 lysosomes/ 
cell. It is encouraging that we obtain from this count a higher num- 
ber than the 60-100 lysosome range we derived in Equation (2), 
in that on the one hand it is plausible to think that not all of the 
lysosomes in a cell may be activated in synchrony, nor may they 
all be involved in NAADP-mediated signaling, and on the other 
hand experience tells us that evolution has built in some redundancy 
of function in order to provide a threshold and also a margin of 
safety for the generation of this type of Ca^+ signals. Moreover, the 
estimated number of junctions is based on a value for [Ca^+Jj^^ 
determined in macrophages, and it is plausible that the total releas- 
able [Ca^^ljy^ in a PASMC may differ from that value and that it may 
also vary over the lysosome population. Lastly, the uncertainty in 



the number of lysosomes/cell emerging from the electron micro- 
graphs as described above may also contribute to this discrepancy. 

L-SR junctional Ca^^ signal 

In the hypothesized model outlined previously ^"'^^ Ca^^ bursts acti- 
vate SR resident RyR3 channels within L-SR junctions and thus 
initiate a propagating Ca^^ wave by CICR. The stochastic simula- 
tions developed here attempt to reproduce the phenomenon of the 
generation of [Ca^^] transients within individual L-SR junctions, 
considering the Ca^^ release kinetic requirements for lysosome- 
resident TPC2 signaling complexes and the rate of Ca^+ capture by 
the SERCA2a localized on the neighbouring SR membrane^^, and 
to link these junctional transients to the observed bursts. 

The thorough quantitative image analysis described in the previ- 
ous section yields critical information for our first modeling phase, 
in which we built a dimensionally accurate virtual lysosome and a 
portion of the SR system, closely apposing the lysosome (Figure 5), 
so as to reproduce a NAADP-triggered Ca^+ signal within the nano- 
space of a representative L-SR junction as faithfully as possible. 

From the available literature we estimated the number of TPC2 and 
SERCA2a likely distributed on the lysosome and SR membranes, 
respectively. We obtained the former number by dividing the mac- 
roscopic whole-lysosome conductance, calculated from the current 
values reported in a recent study on TPC2-mediated Ca^+ current 
in isolated lysosomes^^ by the single channel conductance^\ thus 
estimating that a typical lysosomal membrane may contain » 20 
TPC2. To obtain this value, we used experimental condition data^^ 
to extract values for the ionic potential, E.^^, across the lysosomal 
membrane. We then employed the authors' current- voltage (I-V^ 
data to compute whole-lysosome conductance values (g^) as a 
function of the applied membrane potential, V^, from Ohm's law: 
/ = g^i(y^ - E.^J. We report details on this calculation in the sup- 
plementary material. 

Moreover, we have estimated the density of SERCA2a on the SR 
within L-SR junctions to be equivalent to that previously predicted 
for PM-SR junctions as approximately 6250//im^''^. Other input data 
for the model are the estimates of lysosomal volume and of the 
[Ca^^ljy^, from which we calculated the actual number of ions in 
the lysosome prior to the beginning of Ca^^ release (see previous 
section). This, in turn, was used in the extrapolation of the TPC2 
complex Ca^+ release rate as a function of time, as follows. 

A recent electrophysiological study of the Ca^+ conductance of 
TPC2 signaling complex provides valuable information regarding 
its biophysical properties and, importantly, highlights the fact that, 
for a given relatively low activating concentration of NAADP 
(10 nM, as in Figure 1), the channel open probability appears to 
depend on [Ca^^^^^ (this is likely due to a partial neutralization of 
the electrochemical potential across the lysosomal membrane)"^^. 
We used the Ca^+ conductivity measured in that study and values 
of the lysosome membrane potential^'^ to calculate the maximal 
Ca^^ current, /^^^, as 2.4 x 10^ ions/s. We then applied a weighted 
quadratic fit to the channel's open probability (P^) as a function of 
the reported [Ca^^J^^^ data with constraints that the would tend 
to zero at [Ca^+]j ^ = 80 juM, based on the observation that below 
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Figure 5. A and B, 3D software reproduction of a lysosome closely apposed to a portion of SR, thereby forming an « 20-nm-wide L-SR 
nanojunction; this rendering was inspired by a series of observations from micrographs as in Figure 2 (grey object: SR, blue object: lysosome). 
Included are relevant molecules traversing Brownian motion trajectories produced by the model simulations (see symbol legend below the 
panels). B, enlarged view of the L-SR junctional region, in which we have displayed the volume object (rust-coloured box) used to measure 
the [Ca2+]^j transients like the ones reported in Figure 7. Scale bar: 100 nm. The model geometry and code files are available from the 
corresponding author. 



[Ca^^ljy^ =100 jjM essentially no single channel openings were 
observed^\ Another constraint for the fit was that the curve be 
within the standard deviation value of at the highest reported 
[Ca^^]jy^ = 1 mM. From the quadratic fit, we then obtained our own 
P^-vs-[Ca^+]jy^ table (plotted in Figure 6A) and assumed that at the 
beginning of Ca^^ release, the Ca^^ current would be P x / at 

o o ' o,max max 

the maximal luminal concentration until the luminal concentra- 
tion decreased to the next point in the P^-ys-[CsL^^\^^ relationship. 
At this point, the release rate decreases to a new (lower) P x I 

r ^ ^ ^ o max 

until the luminal concentration reaches the next lower point in the 
table, and so on until = 0 at [C2l^\^ = 80 ^M. The Ca^^ release 
rate as a function of time obtained in this manner is shown in 
Figure 6B. 

To implement an approximated SERCA pump action, we used a 
simplified version of a multi-state kinetic model developed earlier''", 
which was further informed by more recent isoform- specific stud- 
jgg5i,52 particular, since we are interested in simulating the 




0.0001 0.0003 0.0005 

[Ca2-^]iys (M) 

Figure 6. A, open probability of the TPC2-dependent Ca^+ co 
literatures^ . B, Ca^^ release rate calculated as explained in the text. 



SERCA2a Ca^^ uptake only in terms of its influence on the shaping 
of the junctional [Ca^^] transient, we have opted not to implement 
the steps of the multistate model that deal with the Ca^^ unbinding 
from the SERCA on the SR luminal side of the pump. In brief, the 
reactions SERCA2a undergo are: (1) binding/unbinding of the first 
Ca^^; (2) binding/unbinding of the second Ca^^. 

The Ca^^ diffusivity was obtained from studies, in which it was con- 
cluded that given the known kinetics of typical Ca^^ buffers, the 
range of free Ca^+ after it enters a cell can be up to 200 nm. Therefore, 
due to the nano-scale of our system, the trajectories of Ca^^ released 
by TPC2 signaling complexes in the simulations are governed by 
the measured diffusivity of free Ca^^ for a displacement comparable 
to the nanojunction extension (2.23 x 10"^^ mVs;^^'''^). Once Ca^^ 
are buffered they acquire the measured diffusivity of the buffers 
(8.4 X 10-11 m2/s;^^). 

We summarize the quantitative model input data in Table 2. 




time (s) 

nee reproduced from a quadratic fit to the data reported in the 
on the in A. 
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Using MCell as a stochastic particle simulator, we ran a number 
of simulations to represent the NAADP-mediated Ca^+ release 
that is supposed to occur at L-SR junctions in PASM experiments. 
Released Ca^^ is assumed to undergo Brownian motion in the sur- 
rounding space, including the L-SR nanospace. SERCA2a placed 
on the neighbouring SR surface may capture Ca^^ according to our 
approximation of their known multistate model. To determine how 
this regenerated cellular environment can shape a Ca^+ transient we 
"measured" the junctional [Ca^+] by counting the ions within the 
L-SR volume at any given time and dividing the number by that 
volume. The snapshots in Figure 5 are part of the visual output of 
this phase of the work. The data to determine the simulated Ca^^ 
transients in the L-SR junctions were collected in a measuring vol- 
ume placed between the lysosomal and SR membranes in the vir- 
tual L-SR junction (rust-coloured box in Figure 5B). 

Previous experimental findings about VSMC PM-SR junctions 
indicate that disruption of nanojunctions can have profound conse- 
quences on Ca^^ signaling properties of the cells. For example, when 
calyculin A was used to separate superficial SR portions from the 
plasmalemma of the rabbit inferior vena cava, it was observed that 
[Ca^+]. oscillations would cease^^. This was later corroborated by a 
quantitative model of the PM-SR junction's role in the refilling of 
SR Ca^^ during oscillations^'^. In another study by the van Breemen 
laboratory, it was observed that the mitochondria- SR junctions of 
airway SMC displayed a variable average width as a function of 
the state of rest or activation of the cell''^. It makes sense then to 
study the effect of changes in the L-SR junctional geometry on the 
transients generated by our simulations. Therefore, we ran several 
simulations, in which the separation between the lysosomal and SR 
membranes was increased from 10 nm to 100 nm in steps of 10 nm. 
In Figure 7A, we report three sample transients from this set of 
simulations obtained using three different junctional membrane 
separations, as indicated in the inset legend. The value of each of the 



points graphed in this plot is the average value of 100 simulations, 
in each of which the random number generator within MCell is ini- 
tiated with a different seed (see Materials and methods). We report 
representative error bars as 3x the standard error, to convey the 
> 99% confidence interval of the data. As one would expect, the tran- 
sient nanojunctional (NJ) Ca^^ concentration, [Ca^'^Jj^j, decreases in 
magnitude, as the junctional L-SR membrane-membrane separation 
increases, simply because the released Ca^+ has a larger junctional 
volume available over which to spread. However, it is important to 
make a quantitative comparison between this change in [Ca^^J^^j and 
the [Ca^^]. requirements to activate the putative RyR3 population of 
the junctional SR. This can be attained by analyzing the time scale 
of both the recorded and simulated Ca^+ signals. 

Reconciling the temporal scales in simulation and 
experimental results 

In the final step of the development of our model, we analyzed the 
relationship between the time scale of the [Ca^^J^^j transients result- 
ing from our model (Figure 7) and that of the observed Ca^^ signal 
shoulder in Figure 1 . 

Let us note that the typical duration of the simulated transients, 
t . „ is about 0.06 s and recall that these represent A[Ca^+] above 

transient' ^ 

the resting [Ca^^]. On other hand, the build up to the maximum 
value of ArCa^^l ^ ,^ takes about 5 s (we refer to this time as ^ ,^ ; 

-"shoulder ^ shoulder' 

Figure 1). This interval is about two orders of magnitude larger than 
the duration of the simulated individual transients (Figure 7 A). One 
way to reconcile the hypothesis that L-SR junctions are at the base 
of the observed NAADP-induced Ca^^ signals such as the ones in 
Figure 1 and in particular that the signal shoulder emerges from 
lysosomal Ca^^ release at L-SR junctions, is to bring forward the 
assumption that the signal shoulder may be the result of a sequen- 
tial summation effect over many L-SR junctions, each working 
according to an all-or-none mechanism of Ca^^ release, and that the 



Table 2. Quantitative model input data. 



Quantity 


Value 


Source 




300 nM 


[17] and this article 




400-600 AvM 


[44] 


lysosomal average radius 


180 nm 


this article, Figure 3 


(cylindrical) SMC radius 


6 /im 


[46] 


SMC length 


130 ^m 


[46] 


SMC cytosolic volume 


2.4x10-^2 1_ 


see Supplementary material 


TPC2 complexes per lysosome 


20 


[48, 45] and this article 


TPC2 complex Ca^^ release rate 


variable 


this article. Figure 5 


SERCA2a surface density 


6250/A7m2 


[39] 


SERCA2a first Ca^^ binding (unbinding) rate 


5x10«M-is-i (60 s-') 


[50-52] 


SERCA2a second Ca^+ binding (unbinding) rate 


4x^0^W\-'s-' (60 s-^) 


[50-52] 


free Ca^+ diffusivity, D^^^^ 


2.23x10-^° m^/s 


[53, 54] 


Ca2+ buffer diffusivity, D^^^^^^ 


8.4x10-11 m^/s 


[55] 



Summary of input parameters used in various phases of the quantitative model, and references to the origin of their 
values. 
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Figure 7. A, calculated nanojunctional [Ca^+] transient, [Ca^+]^j, "measured" inside the volume of the recreated L-SR nanojunction shown in 
Figure 5. To show the effect of changes in the junctional geometry, we report three transients calculated using different junctional widths of 
20, 50 and 100 nm. B, [Ca^+l^j vs width of junction, concentration values are temporal averages of the transients as in panel A, calculated 
over an interval of 0.065 s (solid circles) and 0.046 s (empty circles); see text for explanation. The shaded area indicates the approximate 
threshold values for CICR at RyRSs^^ 



"firing" of one junction may cause a cascading effect across the set 
of junctions that produces the shoulder. Then the duration of the 
shoulder upstroke can be expressed as: 



^shoulder ^NJ ^ ^transient 



(3) 



where A^^^j is the number of L-SR nanojunctions that yield A[Ca?'^]^^^^^^^^. 
Since we calculated above that Ca^+ release from A^^^j » 60-100 lys- 
osomes would be necessary to produce the observed signal shoulder 
magnitude and have shown in the previous section that such release 
would need to take place at L-SR junctions, we can estimate t . . 

^ ' transient' 

the duration the [Ca^+l^^j transient, by reversing Equation (3) and 
obtaining t^^^^^^^^^ = 0.05-0.08 s. It is noteworthy that this range of 
values is obtained in a manner completely independent of our simu- 
lation results, which yielded a similar range of values. 



To gain quantitative insight into the effects of varying the junctional 
width, we then calculated the temporal average of the [Ca^^]^j over 
^transient (^^i^g ^^c middlc valuc of the range calculated via Equation (3)) 
and graphed it as a function of the junctional width. The result of 
this analysis is reported in Figure 7B (solid blue circles). As we 
have anticipated at the end of the previous section, the decrease in 
magnitude of these data is to be expected, however in this plot we 
also indicated the range of [Ca^^]. values (shaded area) over which 
maximum SR Ca^^ release via RyR3 is reported to occur in skeletal 
muscle^^ This comparison underscores the important constraint 
played by the width of the L-SR junctions and indicates that, unless 
the membrane separation is kept below about 30 nm, it is not possible 
for [Ca^^]j^j to breach the threshold for RyR3 Ca^^ release. Let us also 
point out that the [Ca^^J^^j data in Figure 7B would shift upward, 
toward concentration values that would make the junction more 
prone to promote RyR3 release, if the temporal average were taken 
over a shorter transient time, t' , < t , around the [Ca^+l^,, 

' transient transient -"NJ 

peak. However, in that case Equation (3) indicates that a larger A^^^j 
(than 80, picked as the middle of the 60-100 range) would have 
to contribute to the signal summation that results in a 5-second ^ ,^ . 

shoulder 



Interestingly, this possibility agrees with our lysosome count from 
TEM images (110 lysosomes/cell) and with the argument of natural 
redundancy we contemplated to explain the discrepancy between 
the calculated and observed lysosome number estimates. As an 
exercise we have recalculated the [Ca^^]j^j time-averaged over 
^'transient Obtained uslng 110 lysosomes in Equation (3) (empty pur- 
ple circles in Figure 7B), and this indeed shows that RyR3 Ca^"*^ 
release threshold would be cleared more readily. These observa- 
tions suggest that activation of RyR3s in the L-SR junctions not 
only depends on the concentration of Ca^+ near them, but also on 
their exposure time to this concentration. 



L-SR nanojunction study data 

7 Data Files 

http://dx.doi.org/10.6084/m9.figshare.988699 



Discussion 

Inter-organellar nanospaces 

We have recently introduced the concept of the "panjunctional SR", 
which states that Ca^+ release and uptake at a family of specific 
nanojunctions connected by a continuous but variable SR lumen 
integrates cellular control over multiple functions^. The lysosome- 
sarco/endoplasmic reticulum (L-S/ER) junction is the most current 
junction to be considered in this context and exhibits perhaps the 
highest degree of plasticity of the family of nanojunctions of the 
SR; the mechanism and function of lysosomal Ca^^ signaling is cur- 
rently hotly debated^'^. 

By means of a thorough ultrastructural study in rat pulmonary 
artery smooth muscle, we have observed and characterized L-SR 
nano-j unctions, which had been previously hypothesized on the 
basis of optical measurements of Ca^+ signals and optical immu- 
nocytochemistry"^'^^ Our observations corroborate the previously 
reported finding^^^'^^ that lysosomes in R\SMCs tend to cluster in the 
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perinuclear region, as is evident in our micrographs (e.g., Figure 2A). 
We find that L-SR junctions are on average 15 nm in width (equiva- 
lent to our preliminary reports^'^'" and to recent observations in 
cultured fibroblasts^^) and extend approximately 300 nm in lateral 
dimensions, thereby involving about 15% of the lysosomal mem- 
brane (Table 1). 

Mechanism of NAADP Ca^^ signaling 

In an effort to achieve quantitative understanding of the phenom- 
enon of NAADP-mediated Ca^^ transients and verify the proposal 
that these may be generated in L-SR junctions "^'^^ we focused on 
one of the prominent features of these Ca^^ signals, namely the local- 
ized Ca^+ bursts that precede the propagating Ca^^ wave (Figure 1), 
which we refer to as the signal shoulder (A[C3^^]^^^^^^^^). In the first 
instance, we have estimated the potential contribution of local 
bursts of Ca^+ release from individual lysosomes to the elevation 
of global [Ca^^]. observed during this shoulder in the experimental 
records. From this, and using the dimensions of a typical smooth 
muscle cell, we calculated that 60-100 lysosomes would be 
required to cause an elevation of comparable magnitude to the sig- 
nal shoulder. This is lower than the estimate for the total number of 
lysosomes per cell, 1 10, we obtained from our ultrastructural study, 
but comparable in order of magnitude. We have already mentioned 
above a number of factors in favour of observing a greater number 
of lysosomes/cell than the estimated number required to generate 
the signal shoulder. In addition, several other elements contribute a 
degree of uncertainty to those estimates, such that their discrepancy 
may not be as large. We need to consider that published data 
for the Ca^^ conductance associated with the TPC2 signaling com- 
plex^^ on which we based our TPC2 rate table, show some variabil- 
ity according to the standard deviation bars, which, in turn, implies 
an uncertainty in the interpolated (Figure 5A). Moreover, we 
cannot exclude the possibility that these data reflect a contribution 
from multiple channels (NPJ rather than a purely single channel 
P^. Lastly, the standard deviation of the simulated transient (only 
the standard error is shown in Figure 7A) and the variability in the 
experimental determination of the [Ca^^] sensitivity of RyR may allow 
for some uncertainty in the estimated number of nanoj unctions. 

To take this study further and understand whether the observed 
L-SR junctions could give rise to [Ca^+]. transients of appropri- 
ate magnitude to trigger Ca^^ release from RyR3 channels on the 
junctional SR, we developed a quantitative stochastic model of 
Ca^^ dynamics in the junctional nanospaces. We have previously 
published a simplistic version of this model, which nonetheless 
captured the essential features of the problem and yielded an 
indication that such Ca^^ transients in the L-SR junctions would be 
possible^. However, the simulated transients we obtained displayed 
unphysiological features, such as an abrupt onset and decay. We 
show here that this was largely due to lack of a faithful representa- 
tion of the open probability for the Ca^^ conductance of the TPC2 
signaling complex. We have now combined experimental informa- 
tion on the biophysical properties of conductance and open prob- 
ability"^^ and on the luminal [Ca^^] of the lysosomes^^ to implement 
a more realistic Ca^^ release rate model in the simulations. As a 
consequence, we are able to output a physiologically meaningful 
junctional transient profile (Figure 7A), and observe that the [Ca^+]j^j 
transients generated within our model junctions — in turn, based 



on those observed in our TEM images — reach peaks of about 
20 jjM. While we are still not able to measure these transients in 
individual L-SR junctions (peri-lysosomal Ca^^ probes are only 
recently becoming available^' ^ and so far have not been used in 
smooth muscle cells, vascular or otherwise), it is worth noticing 
that the values we find are comparable to those recently meas- 
ured in so-called Ca^^ hot spots in the mitochondria-ER junctions 
of neonatal ventricular cardiomyocytes (rat culture)^'- and in RBL- 
2H3 and H9c2 cells earlier^''. These results therefore suggest that 
the hypothesis presented for the role of L-SR junctions in cellular 
Ca^^ signaling is certainly plausible. 

Lysosomal trigger zone 

Although our simulations were successful, comparison of the Ca^^ 
transients in Figure 1 and Figure 7 not surprisingly reveals a strik- 
ing difference between the simulated Ca^^ kinetics of a single L-SR 
junction and those of the whole cell. This illustrates important 
aspects related to the concept of lysosome-SR "trigger zone" previ- 
ously introduced by one of us (AME)"^ and recently proposed as a 
facilitating factor in lysosomal-ER signaling in reverse, whereby 
Ca^^ release from ER compartments enables NAADP-mediated 
activation of acidic organelles'-. 

Briefly, this concept suggests that in PSMCs clusters of lysosomes 
and closely apposed SR regions containing sets of RyRs may 
act together to form specialized intracellular compartments that 
transform NAADP- stimulated localized lysosomal Ca^^ release 
into cell- wide Ca^^ signals via RyR- supported CICR. If the firing 
of one junction were sufficient for subsequent initiation of the 
CICR across the entire SR, the experimentally observed threshold 
(Figure 1) would be much lower and the rate of junctional coupling 
by CICR faster. In other words, due to the high [Ca^^]. threshold of 
about 10 for Ca^+ activation of RyR3, CICR engendered by a 
single L-SR junction is likely to die out, unless reinforced by a pro- 
cess of quantal releases by other L-SR junctions within a cluster^'"^. 
Therefore it seems more likely that multiple L-SR junctions work in 
concert in a process characterized by both additive and regenerative 
elements to provide the necessary threshold and margin of safety 
required to ensure the propagation by CICR of the global wave via 
the more widely distributed RyR2 along extra-junctional SR, once 
the Ca^+ release wave escapes an RyR3 -enriched SR region, as pre- 
viously suggested^ ^ 

In this respect, our results also allowed us to establish a time interval 
for the duration of the transient (? . ) that would be compatible 

^ transient^ 

with the hypothesis that sequential summation of all-or-none lyso- 
somal Ca^+ release events from a set of individual L-SR junctions 
is responsible for generating A[Csl^^]^^^^^^^^. Remarkably, this value 
is comparable to the average duration of the simulated transients, 
which was determined on the basis of ultrastructural details and 
completely independently of the summation effect hypothesis. The 
assumption of summation may imply that the role of the NAADP 
as a stimulus is limited to the initial lysosomal Ca^+ release, while 
the recruitment of subsequent L-SR junctions may be governed by 
further lysosomal calcium release, by SR Ca^^ release via Ca^^ 
activated junctional RyR3, and/or propagation and combination 
of calcium signals via inter-junctional clusters of RyR3. There- 
fore, from our model we envision the intriguing possibility of an 
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important regulatory role of lysosomal Ca^+ content by SR Ca^^, in 
such a way that summation of calcium signals at multiple junctional 
complexes may give rise to the shoulder. While only further experi- 
ments can verify this, the quantitative corroboration provided by 
our findings sets on firmer ground the conclusion from earlier stud- 
ies that Ca^^ bursts from lysosomes are indeed responsible for ini- 
tiating the first phase of a cell wide Ca^+ wave via L-SR junctions. 

Disruption of nanojunctions, plasticity and pathology 

Carrying the signal time-scale analysis further and using the cal- 
culated t to take a temporal average of the simulated [Ca^^l^,, 

transient t' o l Jjsjj 

at different junctional widths, we find that above a width of about 
30 nm the transients would be unable to trigger Ca^^ release from 
the RyR3s (Figure 7B). Thus it is possible that heterogeneity and 
plasticity are controlled by a variable width of the junctional nano- 
space. For example, in atrial myocytes it has been proposed that 
NAADP evokes Ca^^ release from an acidic store, which enhances 
general SR Ca^^ release by increasing SR Ca^^ load and activat- 
ing RyR sites^^. This functional variant may be provided by either: 
(1) An increase in junctional distance such that [Ca^^J^^j is insuf- 
ficient to breach the threshold for activation of RyR2, yet sufficient 
to allow for increases in luminal Ca^+ load of the SR via apposing 
SERCA2 clusters; or (2) L-SR junctions in cardiac muscle formed 
between lysosome membranes and closely apposed regions of the 
SR which possess dense SERCA2 clusters and are devoid of RyR2. 
Further ultra- structural studies on cardiac muscle and other cell 
types may therefore provide a greater understanding of how L-SR 
junctions may have evolved to provide for cell- specific modalities 
within the calcium signaling machinery. 

Under metabolic stresses, such as hypoxia, lysosomal pathways 
have been proposed to provide for autophagic glycogen metabolism 
via acid maltase in support of energy supply, and significant lev- 
els of protein breakdown during more prolonged metabolic stress. 
Moreover, although controversial, in some cell types it has also 
been suggested that lysosomes may contribute to the energy sup- 
ply by providing free fatty acids for beta-oxidation by mitochon- 
dria^'^ It may be significant, therefore, that TPC2 gating and thus 
autophagy may be modulated by mTOR '^'^'^\ This may well speak to 
further roles for lysosomal calcium signaling and L-SR junctional 
plasticity during hypoxic pulmonary vasoconstriction and even dur- 
ing the development of hypoxia pulmonary hypertension^^, not least 
because AMP- activated protein kinase has been shown to modulate 
autophagy through the phosphorylation and inhibition of mTOR^l 

The process of autophagy serves not only to regulate programmed 
cell death, but also to recycle organelles, such as mitochondria, 
through a process of degradation involving lysosomal hydrolases'''^. 
For this to occur L-SR junctions would likely be disrupted and thus 
select for local rather than global Ca^^ signals in order to facilitate 
fusion events among lysosomes, endosomes, autophagosomes and 
amphisomes^", which follows from the fact that Ca^^ plays a piv- 
otal role in vesicle trafficking and fusion^ ^ Active non- synchronous 
movements of TPC2-expressing vesicles have been detected in 
live-cell imaging experiments using GFP-tagged proteins-^\ Thus, 
lysosomal Ca^^ signaling may modulate plasticity in the manner 



required at all stages of multiple membrane fusion events that are 
dependent on Ca^^ for the effective formation of the SNARE com- 
plexes^^ In short, transport of proteins between lysosomes, Golgi 
apparatus, and plasma membrane via lysosomes and endosomes 
may be coordinated at different stages by both global Ca^^ signals 
and spatially restricted Ca^^ release from acidic stores in a manner 
that is in some way determined by L-SR junctional integrity. 

Loss of integrity of L-SR junctions may also contribute to disease, 
given that lysosomal Ca^+ release both depletes luminal Ca^+ and 
causes intraluminal alkalinization^-. Therefore, disruption of L-SR 
junctions may modulate the activity of pH-sensitive hydrolytic lys- 
osomal enzymes, such as glucocerebrosidase and acid sphingomy- 
elinase, which exhibit a marked loss of function at pH > 5-^'^' that 
could lead to accumulation of macromolecules such as glucocer- 
ebroside and sphigomylin. As mentioned previously, dysfunction 
of these enzyme systems consequent to L-SR junctional abnormali- 
ties could also contribute to pathologies associated with subclasses 
of lysosome storage disease such as Niemann-Pick disease type 
Cl--'^\ Pompe and Gaucher^^'^^ which may include hepatic por- 
tal-\ or pulmonary hypertension-^-^ dysfunctions in cholesterol 
trafficking^^ and consequent increases in plasma cholesterol levels, 
vascular lesion formation, atherosclerosis/thrombosis and medial 
degradation-^'^^. 

Conclusion 

We have determined that L-SR junctions, about 15 nm in width 
and extending to » 300 nm in lateral dimensions are a regularly 
occurring feature in rat pulmonary arterial smooth muscle cells, in 
which L-SR junctions were first proposed. The present study pro- 
vides a mechanistic basis for the observed NAADP-induced Ca^^ 
signals and strong support for the proposal that these Ca^^ signals 
are generated at L-SR junctions. Even within the variability of our 
recorded values of L-SR junctional widths and extension, our re- 
sults suggest that localized [Ca^^] transients due to junctional Ca^^ 
release can without fail reach values required to breach the thresh- 
old for CICR from junctional RyR3s. Perhaps most significantly, 
disruption of the nanojunctions decreases [Ca^^J^^j below the value 
for CICR via junctional RyR3s. Therefore, consistent with previous 
studies on the PM-SR membrane^'^, we have established that L-SR 
junctions are required to allow such signals to be generated and that 
there is a 30-50 nm limit on junctional width, above which there is 
loss of junctional integrity and inadequate control of ion movements 
within the junctional space. This suggests that the observed L-SR 
junctions in PAMSCs are not only capable of delivering localized 
Ca^+ bursts of the required magnitude, but are also necessary if lys- 
osomes are to fulfill this identified role in Ca^^ signaling. In other 
words, L-SR nanojunctions are a necessary and sufficient condition 
for generating local Ca^^ bursts essential for NAADP-induced Ca^^ 
waves. Most importantly, however, this study demonstrates the 
importance of junctional architecture on the nanoscale to the capac- 
ity for coupling across cytoplasmic nano spaces, tight regulation of 
ion transport and thus signal transduction. In turn, this highlights 
the possibility that alterations in the dimensions and architecture 
of intracellular nanojunctions lead to cell dysfunction and hence 
disease. 
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Supplementary material 



1 Estimating the number of TPC2 complexes per lysosome 

The whole-lysosome conductance was calculated based on the work 
reported by Schieder and collaborators [1], in which current voltage 
relationships of TPC2-mediated Ca^^ currents were determined for 
single lysosomes (see Figure 1 therein [1]). A set of these values, 
as well as the experimental conditions relevant for this calculation, 
are reported in Table SI. From these data, we first calculated the 
electrochemical potentials of (E^ and Ca^^ (E^^ across the lyso- 
somal membrane using the general expression 



find N^p^^ by dividing the calculated conductance, g^, by the sin- 
gle channel conductance, gjp^2- calculation, we picked the 
whole-lysosome conductance value = 316 pS corresponding to 
V = -20 mV, in accordance to the most recent measurements of 

m ' 

lysosomal membrane potential [3]. Hence, our estimate of the num- 
ber of TPC2 signaling complexes in a representative lysosome is: 



^l^=^^2l TPC2 complexes 

^TPC2 1 ^ 



where R is the universal gas constant, T the absolute temperature, 
C^^ is a z-valent cation, F is Faraday's constant. We then obtained 



the whole-lysosome, g^^, 
/ 



from Ohm's law: 



where / is the measured current, is the applied membrane 
voltage, and E.^^ = E^+{PJP^^E^ with PJP^^ the permeability 
ratio. The calculated conductance values are listed in Table SI. 

If we assume that the whole-lysosome conductance under these 
conditions is due to a set of TPC2 complexes working in parallel, 
then = N^p^2 ^tpc2' ^tpc2 number of complexes in that set 
and g is the single TPC2 determined as 15 pS [2]. We can then 



Table S1. Data for whole-lysosome conductance calculation. 



Experimental 
conditions 


V.(mV) 


/(pA) 




i 




60 mM 


100 


15 


265 




2 mM 


60 


5 


301 




70 mM 


20 


-7 


300 




60 mM 


-20 


-20 


316 


T 


296.16 K(23°C) 


-60 


-35 


339 


P IP 


1000 


-100 


-50 


349 


R 


8.3 J/(mol K) 








F 


9.65 X 10^ J/(V mol) 









Data on experimental conditions by Schieder et al. [1], used to calculate 
whole-lysosome TPC2-mediated Ca^^ conductance, g^^, as explained in the 
text. The subscripts "i" and "o" refer to the inside and outside of a lysosome, 
respectively. 



2 Estimating the cytosolic volume of a vascular smooth 
muscle cell 

We model a typical smooth muscle cell as a long and narrow cyl- 
inder 3 jjm in radius and 130 jjm long [4]. Its volume is then V^j^^^ 
= 71(3 /jmY (130 Jjm) » 3.7 x 10"^^ L. To calculate an upper-limit 
estimate of the cell cytosolic volume, V^^^, available to Ca^+ released 
from lysosomes, we need to account for the space occupied by the 
nucleus, the sarcoplasmic reticulum (SR), the mitochondria and 
the set of lysosomes we observe (V 
tively), so that 



V , and V , respec- 

SR' mito' all lys' ^ 



-V -V -V -V ■ 

'^SMC '^nuc '^SR mito 



all lys 



(1) 



The nuclear volume, V , is estimated from a sample of 10 cells in 

' nuc' 

one of our lower resolution micrographs, by measuring the minor 
and major axes of the oval-shaped nuclei, approximating their 
3D shape as a prolate spheroid and calculating the volume of the 
spheroid from the measured values of the axes (Figure 2 and Table 1 
in the main article). We calculate it as 



= 5 X 10"''L 



(2) 



(about 15% of the cell volume). 



To estimate V^.^^, we resort once again to our own micrographs, 
in which we observe of the order of 10 mitochondria per sample 
section per cell. Using a similar procedure to the one employed to 
obtain the approximate number of lysosomes in a pulmonary artery 
smooth muscle cell (described in the main article, section 'Genera- 
tion of the signal's "shoulder" by lysosomal Ca^+ bursts'), and pre- 
tending that mitochondria are spheres 0.8 jum in radius, we arrive at 
an estimate of 43 mitochondria/cell. Thus, 



^™i,„=43x-;r(0.8/^m)^«10-"L 



(3) 
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The volume of the SR in vascular smooth muscle has been estimated 
previously [5] and expressed as a percentage of the cell volume 
deprived of the nuclear and mitochondrial volume. These findings 
suggest that 

y 

-0.75 (4) 

y -y -y 

^ SMC ^ nuc ^ mito 

Which corresponds to 

Ksj,-2 X IQ-^'L (5) 
using the values for V^^,^, V and V . found above. 

^ SMC nuc mi to 

V „ , can be obtained from the value of V, we found in the main 

all lys lys 

article multiplied by the 110 ly so somes/cell we already estimated. 
However, at a value of ~ 3 x 10"^^ L it is negligible when compared 
to the other terms. 

Finally, putting Equation (5), Equation (3) and Equation (2) into 
Equation (1), we find 

F,^=2.4xlO-'^L (6) 

Clearly, the value in (6) is an upper limit, since we haven't accounted 
for all of the organellar volumes in our calculations. However, the 
values above suggest that to a first approximation nucleus, SR, and 
mitochondria are by and large the main contributors to the difference 
between the gross cellular volume and the cytosolic volume. 

3 Line-and-point plot of NAADP-induced Ca^+ fluorescence 
trace 

The plot in the figure below was obtained from the same data, from 
which we produced the Ca^^ trace in Figure 1 in the manuscript. We 
report it again here in line-and-point form to highlight that the sam- 
pling frequency was 0.5 Hz. In a separate supplementary spread- 
sheet file, we also provide the raw data. 
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NAADP mediated two-phase Ca^+ release responses were observed in smooth muscle cells, including 
pulmonary artery smooth muscle cells (PASMCs). Lysosomes and sarcoplasmic reticulum (L-SR) 
junctions had been proposed as the structural bases in facilitating local lysosomal Ca^+ bursts (first 
phase) followed by whole cell Ca^+ release responses from SR (second phase) through a 
calcium-induced calcium release (CICR) mechanism. In their studies reported here, Nicola Fameli etal. 
attempted to define the L-SR junctions in PASMCs with a nano-level resolution using electron microscopy 
(EM), and proceeded to build a quantitative L-SR junctional model. Unfortunately, all these efforts will add 
little to our further understanding of these L-SR junctional structures, if not at all, due to the lack of solid 
experimental EM results in reconstructing L-SR junctions, some dogmatic presumptions in L-SR model 
build up, the inadequacy of validation of the model and a failure to include information on the current 
disputes over the identity of NAADP-associated Ca^+ release channels in lysosomes. Some major 
concerns are as follows: 



1 . Data are not convincing and lack quality. In defining L-SR junctions, the assurance of lysosome 
and SR identities is fundamental. Immune-gold staining has been widely applied in EM studies to 
ascertain different cell organelles as well as to conduct 3D reconstruction of cell structure ^ ^. 
However, such a standard practice was excluded in this study as a 'precaution' to prevent 
compromising lysosome membrane definition. In line with applying empirical judgments for 
identifying organelles to the EM data analysis, the authors extrapolated the opposing and proximal 
organelles adjacent to the nucleus as lysosomes and SR, and deemed them to be worth studying 
(Figure 2). In addition, the authors stated that there were few lysosomes found in the 
subplasmalemmal region, which implied that L-SR junctions were profoundly located around the 
nucleus. However, the fura-2 fluorescence images in Figure 1 as well as in the two citations used 
as basis of this study (10,1 1)^ exclusively demonstrated that the first phase lysosomal Ca^+ 
release initiated from a region around the cell membrane, suggesting that L-SR junctions should 
have cytoplasmic locations near the plasma membrane. It is clear that more reliable methods such 
as immune-gold staining are needed to resolve the contradictions over their own results. 



2. The modeling is based on many assumed parameters and lacks reliable validation. To establish 
the L-SR nanojunctional model, the present study has introduced many assumptions to some 
critical parameters such as the distance between the lysosome and SR, the size of lysosome, and 
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total number of lysosomes; many of which lack proper judgments. In the abstract section, the 
authors rationalized the application of EM in investigating L-SR junctions by stating "L-SR junction 
membrane separation tias been estimated to be < 400 nm and tfius beyond ttie resolution ofligtit 
microscopy". In defining the L-SR nanojunctions, they declared that lysosomes and SR "must be 
separated by a distance of approximately 20 nm and be of the order of a few hundred nm in lateral 
dimensions" by citing PM-SR junctions studies (10,11 )^ ^. However, these two citations originated 
from the light microscope-based studies of L-SR junctions by the same group. Given the resolution 
limitation, these light microscopy results have no way to give an estimation as accurate as 20 nm. It 
is therefore confusing to find that these references fail to back up such claims, in my opinion. In 
reference 39 ^, it has been proposed that the extracellular Ca^+ entered the cell through Na+/ Ca^+ 
exchanges, traversed PM-SR junctions to refill SR via Ca^+-ATPase and sustained cell Ca^+ 
oscillations. For the current L-SR model, this PM-SR modeling may be used as a comparison or 
reference system leading to the conclusion "that L-SR had a similar 20 nm gap measured as that in 
PM-SR junctions." Given that there are no similarities of PM-SR to L-SR junctions, in terms of 
potency of Ca^+ movement due to Ca^+ gradient across PM and lysosomal membrane, 
biophysical property of Na+/ Ca^+ exchangers and TPC2 channels, and the actions of Ca^+ 
-ATPase and RyR3 Ca^+ receptors, such a conclusion is not reliable from the modeling. In 
addition, a major concern is that this model is too simplistic as it excludes many other Ca^+ 
regulatory mechanisms in both lysosomes and SR, which may influence the accuracy of junction 
width or prediction in modeling. 

3. Some assumptions contradict well-established concepts due to a lack of consideration of 
lysosome biology. In explaining the shoulder Ca^+ signal (Figure 1), the authors introduced up to 
four assumptions and calculated that -60-100 lysosomes were needed to generate this shoulder 
peak. They verified these lysosome numbers with their TEM images that "can expect a total of 
about 1 10 lysosomes/ceir. Considering these two calculations however, it is obviously problematic 
because almost all lysosomes (100 out of 1 10 each cell) will need to functionally contribute to the 
L-SR junctional structure formation. In addition, TEM measurements showed the lysosome 
dimensions (Table 1) to be major axis of 398 ± 91 nm and minor axis of 325 ± 65 nm. Given the 
L-SR junction model features "junctional membrane separation and extension are about 15 nm and 
300 nm", it could be postulated that all the lysosomes that participated in L-SR structure formation 
also have to align well with the SR in order to produce a 300 nm junctional extension. Therefore, 
this L-SR model actually describes a scenario whereby in PASMCs the majority of lysosomes act 
as the components of L-SR junctions and uniformly have impeccable arrangements to oppose the 
SR at 20 nm apart. However, this idealistic model largely contradicts the well-established concepts 
that lysosomes are active organelles that dynamically exert their different functions by their 
movements and other activities. 

4. TPOs as lysosomal Ca^+ release channels are questionable. In quantifying the L-SR model, the 
authors arbitrarily assumed TPC2 as the only lysosomal NADDP-activated Ca^^ channel by totally 
disregarding recent findings that TPC2 is not a lysosomal NAADP-associated Ca^+ channel, but a 
Na+ channel^, . They also ignored the contribution of other lysosomal channels or transporters 
that may contribute to such L-SR junction activity. This prejudice in channel selection inevitably 
invites doubts on the reliability of this proposed L-SR model. In addition, the authors collectively 
adopted electrophysiological data from other research groups without taking into consideration any 
incompatibilities with the model. For example, in obtaining the total number of TPC2 channels on a 
lysosome, the authors "...dividing the macroscopic whole lysosome conductance, calculated from 
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the current values reported in a recent study on TPC2-mediated Ca^"" current in isolated 
iysosomes (48) ^, by ttie single channel conductances(45) ^, thus estimating that a typical 
lysosomal membrane may contain --20 TPC2". Unfortunately, the Iysosomes used in the whole 
lysosome conductance recording were not normal Iysosomes, but Iysosomes manipulated by 
overexpression of TPC2 and were undergoing size enlargement. With these changes, the modified 
Iysosomes would have a much higher TPC2 conductance compared with normal Iysosomes with a 
natural biogenesis. It is important to use a TPC2 conductance and other parameters of channels 
from Iysosomes of normal cells without an artificial transgene for such modeling. Otherwise, the 
modeling prediction will be inaccurate and the relevance of such a model will be vague. 

In conclusion; aside from the issue that the current L-SR nanojunction model fails to contribute any new 
insights to our understanding of the L-SR architecture, my concerns over the reliability of data and 
modelling parameters are rather serious - this work is therefore 'Not Approved' in its current version. 
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I have read this submission. I believe that I have an appropriate level of expertise to state that I 
do not consider it to be of an acceptable scientific standard, for reasons outlined above. 
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This is certainly an interesting study. The characterisation of lysosomal-SR junctions is novel and 
conceptually important. The data look convincing but I would recommend the authors to further expand 
the structural part of the study. It would perhaps be useful to further verify that the organelles participating 
in the formation of the contacts are indeed lysosomes and that the junctions are not simply the result of 
accidental collisions of the two organelles. Identifying/visualising tethers between the organelles (by 
analogy with tethers between mitochondria and the ER'') would be particularly advantageous. 

Modelling part of the paper is also useful. In this respect it would be valuable to clarify the Ca^+ buffering 
capacity of lysosomes (and its possible effect on the outcomes of the model). One could also perhaps 
consider measuring this important parameter using isolated lysosomes. 
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This is an interesting and thought provoking study which simulates NAADP-dependent Ca^+ transients; 
these are produced by release of Ca^+ from lysosome stores that induce Ca^+-induced Ca^+ release from 
SR stores via stimulation of RyR3. It is increasingly apparent that answers to many issues concerning 
cellular Ca^+ signalling will be found by recording local Ca^+ events - for instance the recently identified 
Ca^+ sparklets representing Ca^+ influx through opening of plasma membrane ion channels. 

Simulation studies are always difficult to review and evaluate as they must make many assumptions. 
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However, this work is considered and may be a useful model for how other local Ca^^ events might be 
modeled between discrete organelles. Moreover, unlike many simulation manuscripts, I find this one clear 
and well-written. 



Do the authors think that the use of proximity ligation assay techniques might be appropriate to measure 
local interactions between TPC2 and RyR3? Also, how are lysosomes refilled with Ca^+ following 
NAADP/TPC2 activation - a role for a lysosome store-operated Ca^+ influx? 
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